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THE  EFFECT  OF  NOISE  ANISOTROPY  ON 
DETECTABILITY  IN  AN  OPTIMUM  ARRAY  PROCESSOR 

Franz.  B.  Tuteur 

Department  ot  Engineering  and  Applied  Science 
Yale  University 
New  Haver.,  Connecticut  06520 


ABSTRACT 

The  effect  of  localized  noise  sources  on  the  perform¬ 
ance  of  the  optimum,  i.e.,  likelihood  ratio  detector,  is  in¬ 
vestigated. 

Expressions  were  obtained  for  the  performance  loss 
of  optimum  detectors  with  noise  which  consists  of  (1)  an 
isotropic  part,  and  (2)  a  component  from  multiple  point 
source  b. 


INTRODUCTION 

The  object  cf  this  paper  is  to  investigate  the  effect  of  localized  noise  sources  on  the  per¬ 
formance  of  the  optimum  (likelihood- ratio)  detector. 

In  previous  work1  expressions  for  the  performance  loss  of  likelihood- ratio  detectors  when 
the  noise  consisted  of  an  isotropic  part  and  a  component  from  a  single  point  source  were  ob¬ 
tained.  In  the  present  analysis  these  results  are  extended  to  the  case  of  more  than  one  point 
source  in  an  attempt  to  also  get  some  estimate  of  the  performance  loss  caused  by  anisotropy 
sources  that  are  not  strongly  localized.  Such  sources  can,  presumably,  be  represented  by  a 
large  number  of  closely  spaced  point  sources. 


NOMENCLATURE 

The  notation  used  is  similar  to  that  used  by  Edelblute,  et  al.2  The  detector  is  assumed  to 
be  a  directional  array  consisting  of  M  hydrophones,  and  the  received  signal  at  the  ith  hydro¬ 
phone  is  xjft).  Then  if  the  spectrum  of  x;(t)  is  limited  to  frequencies  below  W  Hz,  and  the 
*(t)  are  observed  over  an  interval,  T,  such  that  wt  »  I,  x;(t )  can  be  expanded  in  a  Fourier 
Series: 


1Peter  M.  Schultheiss,  "Passive  Detection  of  a  Sonar  Target  in  a  Background  of  Ambient  Noise 
and  Interference  from  a  Second  Target,"  Yale  University  Progress  Report  No.  17,  submitted 
to  General  Dynamics  Corporation,  Electric  Boat  Division,  September  1964. 

This  report  is  included  in  "Processing  of  Data  from  oonar  Systems,  Volume  111,"  Kanefsky, 
Levesque,  Schultheiss,  and  Tuteur,  General  Dynamics  Corporation,  Electric  Boat  Division  Re¬ 
port  U417-65-033  (Aug.  23.  1965). 

^David  J.  Edelblute,  Joanne  M.  Fisk,  and  Gerald  L,  Kinnison,  "Criteria  for  Optimum- Signal - 
Detection  Theory  for  Arrays,"  J.  Acoust.  Soc.  Am.,  41,  199-205  (Jan.  1967). 
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*l(t) 


=  xi(n)  * 


j  Jnnt/T 


(i) 


where  the  x4(n)  are  complex  Fourier  coefficients  satisfying  x;(-n)  =  x|(n)  and  where  the 
asterisk  stands  for  complex  conjugate.  All  the  available  information  about  the  signals  received 
by  the  entire  array  is  therefore  contained  in  the  set  of  vectors 


X(n)  = 


x,(n) 


(2) 


It  Is  assumed  that  X(n)  and  X(m)  are  statistically  independent  for  n  /  ±m.  Suppose  that  the  sig¬ 
nal  xt(t)  received  at  the  ith  hydrophone  consists  of  signal  and  noise;  then  the  signal  component 
is  gtven  by 


■■<»  •  L 


yi(n> 


i  I»n*/T 


(3) 


so  that  the  signal  component  at  all  hydrophones  is  represented  by 

r 


Y(")  = 


(4) 


Ly»<n>J 

Here  again  Y(n)  are  assumed  to  be  independent  from  Y(m)  for  n  t  tm.  Also,  the  signal  is  as¬ 
sumed  to  be  independent  from  the  noise.  The  normalized  noise  covariance  matrix  is  defined  by 


5<n)  =  nJtT)  <-’<n)  5T(n)>N  •  (5) 

where  the  superscript  T  refers  to  matrix  transposition  and  the  symlxil  <  >N  means  ensemble 
average  subject  to  the  noise-only  hypothesis.  N(n)  is  the  average  noise  power  at  frequency 
2im/T  radians  per  second. 

The  normalized  signal  covariance  matrix  is 


-(n)  =  s<!o  <X*(n>  -T(n)>s  ’  (6) 

where  S(n)  is  the  average  signal  power  at  frequency  tm/t  radians  per  second.  If  the  signal  is 
a  plane  wave,  the  elements  y4(n)  of  Y(n)  are  all  delayed  replicas  of  each  other;  thus 


y4(n)  =  C4*(n)  e  T  ,  (1) 

where  s(n)  is  the  nth  Fourier  coefficient  of  the  signal  wave  form;  the  c.  are  weighting  factors 
to  take  into  account  that  the  signal  strength  or  gain  at  different  hydrophones  may  be  different, 
and  Tj  is  the  delay  at  the  ith  hydrophone.  The  c4's  are  conveniently  defined  in  such  a  way 
that 


<s»(n)  s(n)>  =  1 


(8) 


for  all  n.  Hence 
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and  therefore 


Y(n)  =  s(n) 


s(n)  V„(n) 


S'">  ■  «<»»  !.’(•) 


P(n)  =~Y;(n)  V0T(n). 
S(n) 


(9) 


UO) 

(U) 


P(n)  is  seen  to  be  of  rank  1 .  Because  of  the  independence  of  signal  and  noise  the  covariance 
matrix  of  signal  and  noise  together  is  N(n)  Q(n)  +  S(n)  P(n). 

The  detection  performance  of  the  optimal  processor  is  defined  in  terms  of  the  standard 
detection  index.3 


d  = 


(12) 


where  m  ,  and  m0  are  the  mean  values  of  the  output  signal  when  signal  is  present,  and  absent, 
respectively,  and  where  is  the  standard  deviation  of  the  output  under  the  condition  that  sig¬ 
nal  is  absent. 

By  means  of  a  trivial  extension  of  the  result  of  Edelblute  et  al.,2  it  can  be  shown  that 


d  = 


£  K(n)  S(n)  C0J(n) 


-y/^  K2(n)  N 3(n 


)  G0J  (n) 


(13) 


where  G0(n)  is  the  maximum  value  of  the  array  gain  at  frequency  2nn/T,  given  by 


G0(n)  -  V0T(n)  Q‘l(n)  Vj(n)  , 


(14) 


and  where 


-  S(n)/N»(n) 

1  +  S(n)  0#(n)/N(n)  ' 

Equation  (13)  can  be  simplified  somewhat  by  using  a  small  signal  approximation:  if 


(15) 


3F.  Bryn,  "Optimum  Signal  Processing  of  Three-Dimensional  Arrays  Operating  on  Gaussian 
Signals  and  Noise,"  J.  Acoust.  Soc-  Am.,  34,  289  (Mar.  1962). 
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S(n)  G0(n)  «  N(n)  . 


then 


and 


K(n) 


S(n) 

NJ(n) 


S2(»)  C„(n)2 
NJ(n) 


(10) 


(17) 


(18) 


Equation  (18)  will  be  employed  in  the  sequel,  under  the  assumption  that  inequality  (16)  usually 
holds. 


EFFECT  OF  DIRECTIONAL  INTERFERENCE 

Suppose  that  the  noise  component  of  x.(t)  consists  of  two  parts,  an  Isotropic  part  and  an 
interference  part.  It  is  assumed  that  the  interference  is  generated  by  R  point  sources.  The 
rth  point  source  is  located  at  an  azimuth  angle  er,  and  its  spectral  density  is  !,(&>);  hence  the 
interference  power  from  the  rth  interference  source  at  the  frequency  un  is  ir(n).  The  desired 
target  Is  at  the  azimuth  angle  6t  -  o,  and  it  is  assumed  that  the  array  is  steered  in  the  target 
direction.  The  isotropic  noise  power  at  the  frequency  o.n  is  N0(n).  The  isotropic  noise  compo¬ 
nent,  the  Interference  sources,  and  the  target  signal  are  all  assumed  to  be  mutually  independent 
Gaussian  processes  with  zero  mean.  Then  the  total  noise  power  density  is  given  by 

R 

N(n)  =  N0(n)  +  £  U")  .  (19) 

■  *  l 


and  the  normalized  noise  covariance  matrix  has  the  form 


Q(n> 


N0(n) 

N(n) 


V”’  1  r  (n  )  T 

9o(">  'L  NOO  X.  (n)  ' 


(20) 


where  Q#(n)  is  the  normalized  covariance  matrix  of  the  isotropic  noise  component  and  where 
each  element  of  the  summand  results  from  one  of  the  interference  point  sources.  By  direct 
analogy  to  Eqs.  (9),  (10),  and  (11) 


<  o 


Yr(n)  = 


(21) 


< ') 
2rnrH 


where  T(1r)  is  the  delay  of  the  plane  wave  from  the  rth  interference  source  at  the  ith  hydro¬ 
phone. 
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The  matrix  Q(n)  can  be  Inverted  by  using  the  following  matrix  identity:  il  A  is  a  non¬ 
singular  matrix  of  dimension  M  and  B  is  a  matrix  of  M  rows  and  R  columns,  then 

(A  1  5*  5T)‘  '  =  A'  1  -  A" 1  B*  (I  t  Bt  A'  1 U* )*  ‘  BT  A*  1  (22) 

This  identity  is  easily  proved  by  multiplication. 

In  the  present  application  let 


N.(n) 


(23) 


and 


B* 


vWf v-'(n) 


1 


N(n)  *">J 


(24) 


Also,  to  simplify  the  notation  let 


Kr  =  Kr(n)  = 


I,(o) 

N(o) 


for  r  =  1 


(25) 


and 


G,t  =  «„(")  =  Y,T(n)  Oo  *(n)  Y>)  . 


(26) 


Note  that  Gr,(n)  has  the  general  form  of  an  array  gain  (see  Eq.  (19));  we  can  call  it  a  "cross 
array  gain."  It  is  clear  that 


G*,(n)  =  G„(n>  .  (27) 

In  terms  of  this  notation  the  matrix  I  +  §T  A* 1 B*  of  Eq.  (22)  becomes 


1  ♦  K.G,, 

|  Kj  Gjj  *  * 

•  v’Ki  Ka  gir 

VkTkT  i 

1  ♦  KaG2a  .  . 

•  v  Kj  GjR 

"  I  +  G 

(28) 

g*ir  . 

.  1  ♦  K„  gbr 

by  an  obvious  definition  of  the  "cross-array-gain"  matrix  G.  Note  that  G  is  square,  of  dimen¬ 
sionality  R,  and  Hermitian. 

By  defining  a  vector  of  cross- array  gains  g  such  that 

gT  =  Go  i  V k"a"  co  i  ■  •  ■  vTg  Gor  | ,  (29) 

and  detection  index  d,  as  given  in  Eq.  (18)  can  be  written  in  the  form 
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K(»)  1  ,Yo‘0'>  0‘ '('»)  v;(n) 


-  Goo  -  KTlI  +  g]*,R*.  (31) 

It  Is  clear  that  G00  1b  the  optimum  array  gain  with  Isotropic  noise  and  that  gTll  +  c)' '  k* 
represents  the  effect  of  the  interference.  In  general,  the  evaluation  of  the  interference  term 
In  specific  instances  is  difficult  for  two  reasons: 

1.  The  cross-array  gains,  G,  j(  arc  quadratic  forms  Involving  the  inverted  Q„(n)  matrix. 

2.  Even  If  the  G,  j  are  known  the  R »  R  matrix  [I  ♦  g)  must  be  inverted.  Thus  it  is  neces¬ 
sary  either  to  solve  Eq.  (31)  by  computer  or  to  make  approximations  permitting  an  analytic 
result. 

The  standard  simplification  that  has  been  used  in  most  previous  analyses  and  which  elimi¬ 
nates  part  of  the  difficulty  involved  in  evaluating  the  G, ,  Is  to  assume  that  there  Is  no  correla¬ 
tion  between  different  hydrophones  due  to  the  isotropic  noise  component.  The  effect  of  correla¬ 
tion  has  been  considered  before  and  is  not  expected  to  alter  the  results  obtained  here  in  any 
significant  way.  with  this  assumption  Q0(n)  -  I ,  and 


«.*(">  =  tlu(n)  -  =  Ggg(n)  =  M 

JL,  .  .  (r)  (•) 

«..(»)■  r  ”•  ’• 

TT| 


(32) 

(33) 


A  further  small  simplification  results  from  the  assumption  that  the  array  Is  steered  on  target; 
this  implies  that  t|#>  =  o,  and  therefore 


G0,(n> 


(34) 


Equation  (31)  is  now  explicitly  evaluated  for  a  number  of  simple  special  cases. 


SINGLE  POINT  INTERFERENCE 

If  there  is  only  a  single  interference,  Eq.  (31)  becomes 

K,lG0,|J 

Go°  -  1  "k7g~  ’ 

and  with  the  simplification  Q0(n)  =  i;  this  becomes 

K.iCo.l* 


F(n)  --  M 


1  +k,M  * 


by  Eq.  (34) 


(35) 


(38) 
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If  the  two  interference  sources  are  widely  separated  in  angle  from  each  other  and  from  the 
target,  then  rt  *>  and  r<  J>  differs  substantially  for  all  k  and  are  not  close  to  zero.  Under  these 
conditions  the  coefficient  of  k,k 2  in  both  the  numerator  and  denominator  of  Eq.  (40)  increases 
with  M  while  the  other  terms  increase  with  NJ.  Hence,  for  large  M  these  coefficients  become 
negligible  with  the  result; 


F(n) 


M  - 


X,|G0llJ  Kj|G0  j| : 


and 


F(n)  %  M 


1  t  K,H  1  i  K2H 


KtM 


1  t  K.M 


(42) 


(43) 


where  the  second  approximation  Involves  neglect  of  the  oscillating  terms  in  |  G0 ,  | 2  and  |  G0  2 1 2 . 
Thus  the  effect  of  interferences  is  seen  to  be  additive  under  these  conditions.  For  small 
interference- to-^imbient-noise  ratios,  where  K,M  and  K2M  are  very  much  less  than  unity,  F(n) 
is  reduced  roughly  by  (K,  +  K2)M ;  while  for  very  large  interference-to-amblent-noise  ratio,  the 
reduction  is  no  greater  than  2.  Thus  for  small  interference-to-noise  ratio,  the  detection  index 
d  decreases  roughly  with  the  first  power  of  interference-to-noise  ratio  (see  Eq.  (29)),  but  the 
maximum  effect  is  no  greater  than  the  loss  of  two  hydrophones. 

Suppose  next  that  the  two  interference  sources  are  sufficiently  close  together  so  that  for  all 
frequences  of  interest,  and  for  all  i. 


then 


-n|r!‘>  -  r<2>|  *0, 


K.I*. 


(44) 


and 

Then 

F(n)  %  G00  — 


®e  (Gj,  G,jGj  j)  ~  M|G0  ,| 1  . 

K,{1  +  KjM)|G01  | 1  +  K,<1  +  K,II)|G01|  2  -  2MK tK 2 1 G0  1 1  2 
1  +  (K,  +  K,)M 


(45) 


(K,  ♦  K;)|G0  ,|  1 

Gqo  ,  .  (46) 


00  1  +  (K,  +  Kj)H 


Thus,  the  result  converges  to  the  case  of  a  single  plane-wave  interference  of  strength  K,  *K, 
in  this  case. 


For  a  linear  array  of  M  elements  spaced  d  feet  apart  we  can  set 

|t(.  l)  -  t(,2)|  =  j  1 11/2  -  i|  J  I  sin  0,  -  sin  0,||  , 


(47) 
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where  the  delay  at  the  central  array  element  Is  arbitrarily  given  the  value  zero.  For  small 

ei "  ei> 

|r<‘)  -  rj*)|  =  2  |  M/2  -  i|  ^  sin  |  j  cos  *  I  M/2  -  *1  f  cos  .  (48) 

where  =  (9,  +  <?3)/2  is  the  angle  hallway  between  the  two  interferences.  Then,  since  wn 
iS  2trW 

“n|TI,>  '  Ti2>l„..  =  nWM  ?  1^1  *0ll  eos  6m  (49) 

and  the  two  interferences  are  close  enough  together  so  that  Eq.  (46)  holds  if 

7tWM  £  |0,-02|  cos  en  «  1 


As  an  example  let  II  =  5000  Hz,  d  =  2  ft,  c  =  5000/ft,  and  M  =  10,  then  if  |0,  -  8} \  «  o. 016/cos  da 
radians  the  two  interference  points  have  the  same  effect  as  a  single  one  with  a  higher  power 
level  and  therefore  the  maximum  detectability  loss  can  be  no  greater  than  a  single  hydrophone 
as  shown  by  Eq.  (39). 

Note  that  for  6n  =  v/2  the  two  interferences  are  located  symmetrically  to  the  end-fire  axis 
of  the  array;  therefore  their  effect  is  always  that  of  a  single  interference.  This,  however,  is 
due  to  the  symmetry  of  the  linear  array  and  does  not  hold  in  other  cases. 

Note  further  that  Eq.  (50)  is  a  rather  conservative  limit  since  neither  the  effect  of  integra¬ 
tion  over  frequency  or  over  hydrophone  spacing  has  been  considered.  Depending  on  the  exact 
form  of  the  power  spectrum  these  integrations  should  result  in  increasing  the  value  of  |0,  -  02\ 
by  a  factor  of  4  or  5  over  that  given  in  Eq.  (50). 

If  the  interference-to-noise  ratio  is  small  enough  so  that  (K,  *  K3)  M  <<  1,  then  Eq.  (46)  and 
Eq.  (43)  are  approximately  the  same;  thus  under  this  condition  the  effect  of  cwo  interferences 
on  the  detectability  is  proportional  to  the  interference  power,  and  independent  of  the  spacing  of 
the  two  interference  sources  from  each  other.  Note,  however,  that  approximating  Ig01!  2  and 
I G0  2 1 2  by  M  still  implies  that  the  interference  direction  is  substantially  different  from  the  tar¬ 
get  direction. 


MORE  THAN  TWO  INTERFERENCES 

Extension  of  the  results  obtained  so  far  to  more  interferences  is  difficult  and  involves  fur¬ 
ther  approximations.  Consider  first  the  large  Interference-to-noise  ratio  case,  with  all  inter¬ 
ferences  widely  separated.  We  assume  as  before  that  Q0  =  I  and  that  therefore  Gtr  -  M  for 
r  =  0. 1.2. . .  ,R.  It  can  be  seen  from  Eq.  (41)  that  the  off-diagonal  elements  of  G  are  of  order 
\4F.  It  can  be  shown  in  general  (see  Appendix)  that  for  the  purpose  of  approximate  inversion  an 
n-dimensional  matrix  whose  diagonal  elements  are  of  order  k  relative  to  the  off-diagonal  ele¬ 
ments  can  be  approximated  by  a  diagonal  matrix  if  k  »  n.  Thus  if  i/H  R  the  off-diagonal 
elements  of  the  matrix  [I  +  G]  can  be  neglected  in  forming  the  inverse,  with  the  result  that 
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F(n)  %  G00 


Y  KrlGr0|2 

Z_i  1  +  Kr  G 

r c  1  r  r  r 


(51) 


r~i  KM 

"'Litm 

»*  I 


+  IZ]  H  oos  w" 

»  3  1  k»  i  ♦  1 


(52) 


where  T<r>  is  the  interference  delay  from  the  rth  interference  source  at  the  ith  hydrophone. 
Equation  (51)  is  the  direct  extension  of  Eq.  (43)  and  indicates  that  for  widely  spaced  point  inter¬ 
ference  sources  the  detectability  loss  is  approximately  equivalent  to  the  loss  of  one  hydrophone 
per  interference  source.  The  approximation  is  good  only  for  R  « 


EFFECT  OF  DISTRIBUTED  INTERFERENCE  SOURCE 

A  distributed  interference  source  can  be  represented  by  a  large  number  of  closely  spaced 
point  sources.  Suppose  that  the  interference  source  has  a  spectral  density  l(n>  and  that  the 
interference  power  is  uniform  for  angles  inside  the  interval  <  6  <  92  and  zero  outside. 

Then  the  interference  can  be  represented  by  R  points  of  spectral  density  i(n)/S  equally  spaced 
in  the  interval,  where  r  is  a  large  number.  Initially  it  will  be  assumed  that  the  interference- 
to-ambient-nolse  ratio  is  small.  Although  the  result  obtained  under  this  assumption  is  some¬ 
what  academic  (since  the  interference  effect  is  very  small  in  any  case)  it  is  possible  to  obtain 
an  analytic  result  which  is  probably  applicable  with  some  mofidications  to  larger  interference- 
to- ambient- noise  ratio  as  well.  Under  this  assumption,  the  elements  of  the  matrix  G  are  all 
very  small,  and  it  is  approximately  true  that 

I  +  G  %  I  .  (53) 

Then  the  matrix  inversion  is,  of  course,  trivial.  The  precise  condition  for  Eq.  (53)  to  be  a  good 
approximation  may  be  deduced  from  Ref.  (4);  a  simple  sufficient  condition  is  that 

R 

I0"1  «  1  for  ail  r  -  1,  .  . . ,  R  .  (54) 


In  the  present  discussion  Kr  =  K,/R  for  all  r  =  1 . R,  where  Kx  =  I(n)/N„(n)  is  the  total 

interference-to-ambient  noise  ratio.  A  conservative  upper  bound  on  Kj  such  that  Eqs.  (53)  and 
(54)  are  good  approximations  is  obtained  by  letting  |cr,|  -  u  for  all  r.s  (see  Eq.  (49)).  Hence, 
if 


KZM  «  1  , 

Eq.  (53)  is  a  good  approximation,  and  under  these  conditions  Eq.  (31)  becomes 

F(n>  =  Gao  ~  ET1* 


(55) 


(56) 


and  by  use  of  Eq.  (41)  this  becomes: 


^B.  Freeman,  Principles  and  Techniques  of  Applied  Mathematics  (John  Wiley  and  Sons,  New 
York,  1957),  p.  34. 
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F(n) 


°oo  -  x,m  -  v  £  L tos  “’n  (T*r) 

i=l  k=  i  ♦  1  f*i 


(57) 


We  assume  now  that  the  azimuth  angle  subtended  by  the  interference  is  small  enough  so 
that  the  do  not  differ  very  greatly  as  r  goes  from  l  to  R.  Then,  it  is  possible  to  expand 
r(r>  in  a  Taylor  series  in  r  as  follows: 


r<r>  -  r(m>  +  (r-in)  ,  (58) 

where  m  =  R/2  is  used  as  the  point  about  which  the  expansion  is  performed;  T<m>  is  effectively 
the  mean  delay  of  the  interference  wavefront.  1 

As  R  is  allowed  to  go  to  infinity  the  summation  in  r  can  be  converted  into  an  integral. 
Hence  Eq.  (57)  becomes,  after  some  reduction 


where,  in  going  from  Eq.  (59)  to  (60)  we  have  used  the  fact  that 

f  (Atj  -  Ark)  =  m(ATj  -  Ark)  =  T<m>  -  r<n>  .  (61) 

As  before,  the  term  representing  the  loss  of  detectability  is  the  bracketed  term  in  Eq.  (60). 
Except  for  the  sin  x/x  term  the  form  of  the  double  summation  is  the  same  as  that  which  would 
be  obtained  for  a  single  interference,  (Eq.  54)  with  mean  delay  at  the  ith  hydrophone, 
under  the  condition  K,H  «  1.  In  fact,  the  argument  that  the  summation  of  the  oscillating  terms 
tends  to  become  negligible  applies  here  with  ever  greater  force,  because  of  the  sin  x/x  term. 
One  can  conclude,  therefore,  that  for  interference- to- ambient  noise  ratio  small  enough  to  sat¬ 
isfy  Eq.  (55),  and  if  the  angle  subtended  by  the  interference  is  relatively  small,  a  distributed 
interference  source  affects  the  performance  in  essentially  the  same  way  as  a  single  point 
interference. 

In  order  to  obtain  an  rstimate  of  the  magnitude  of  azimuth  angle  that  can  be  considered 
"small,"  consider  a  linear  array  with  M  hydrophones  spaced  d  feet  apart.  For  such  an  array 

r!r)  =  i  -  sin  er  ,  (62) 

1  C  r 


where  8  r  is  the  azimuth  angle  of  the  rth  interference  point.  Assume  that  the  interference 
power  is  uniform  over  the  range  8X  <  8  <e2  and  is  zero  outside  this  range.  The  center  of  the 
interference  is  at  the  angle 


«.=  {<*,  +  *,)• 

Then,  by  analogy  with  Eq.  (58)  we  expand  sin  0r  about  8m;  i.e., 


sin  8r  *  sin  8m  +  (0,  -  0n)  cos  8m  . 


(63) 


(64) 
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All  the  other  steps  leading  to  Eq.  (60)  can  then  be  performed  In  exactly  the  same  way,  with  the 
summation  over  R  replaced  by  an  Integration  over  ot .  The  final  result  an  be  put  Into  the  form 


F(n)  -  G00  K| 


H  +  2 


-  .  fkwnd  . 

ST^  (  ko)nd  \  Sin  -J-  COS  0m  I  - r —  I 

/  (M-  k)  cos  I  -  sin  8  I - ■  ^  —  - t-t3 —  . 

/  ,  \  c  m>  ,  /<VM 

kTi  —  cos  <?„  ^ - —  ) 


(65) 


Except  for  the  sin  x/x  term  in  the  summation,  this  is  again  the  expression  that  one  would  have 
obtained  for  a  single  point  interference  location  at  the  angle  0m.  It  is  clear  that  the  accuracy 
of  this  expression  depends  on  the  accuracy  of  Eq.  (64)  which  in  turn  is  a  fairly  good  approxima¬ 
tion  for  &,  -  9 ,  less  than  about  1  radian.  Thus  we  conclude  that  an  Interference  source  spread 
over  no  more  than  cne  radian  affects  the  detectability  essentially  like  a  single  point  interfer¬ 
ence  provided  that  KZM  «  1. 


Since  the  effect  of  interference  for  K,M  «  l  is  very  small,  the  result  )ust  obtained  is 
somewhat  academic  and  it  would  be  desirable  to  extend  it  somehow  to  the  case  of  KjM  »  1. 
Unfortunately  this  is  quite  difficult;  in  fact,  the  only  simple  result  that  has  been  obtained  is  an 
extension  of  Eq.  (46)  to  more  than  two  interference  sources.  As  in  the  case  of  two  interferences, 
it  is  assumed  that  the  interference  points  are  close  enough  together  so  that  for  all  frequencies 
of  interest 


max  wn  -  -r**1]  0  .  (66) 

The  extension  to  R  interference  points  is  then  quite  straightforward,  and  the  result  is  that  the 
detectability  loss  is  again  equivalent  to  that  of  a  single  Interference  source  of  strength  K, .  For 
a  linear  array  having  M  equally  spaced  hydrophones  the  maximum  value  of  0,  -  0,  for  which 
this  result  holds  is  given  by  Eq.  (50). 


COMPUTATIONAL  RESULTS 

Since  it  has  not  been  possible  to  obtain  meaningful  analytic  results  for  cajes  in  which  the 
approximations  made  in  the  above  work  are  not  applicable,  F(n)  has  been  evaluated  on  a  digital 
computer  for  a  number  of  different  array  and  interference  patterns,  and  for  specific  frequen¬ 
cies.  The  results  of  some  of  these  computations  are  presented  in  Figs.  1  through  6.  In  all 
computations  it  is  assumed  that  the  array  is  steered  on  target  at  an  angle  8  -  o  and  that  an 
interference  exists  at  some  angle  0,.  The  curves  are  then  plots  of  F(n)  as  0,  is  varied.  Thus, 
if  0,  is  near  o  the  interference  is  near  the  target  in  azimuth,  and  F(n)  is  small.  Also,  the  as¬ 
sumption  that  correlation  of  ambient  noise  waveforms  between  different  hydrophones  is  zero 
has  not  been  used;  instead  the  exact  form  of  the  Q„(n)  matrix  as  given  by  Bryn*  was  used.  As 
a  result  F(n)  *  M  in  the  absence  of  interference  as  would  be  Inferred  from  equations  such  as 
(36),  (43),  and  (52).  In  fact,  F(n)  <  M  in  all  cases;  however,  this  is  a  coincidence.  It  is  possible 
for  F(n)  >  M  as  is  shown  by  Bryn.3  In  all  cases  die  interference-to- ambient  noise  ratio  is 
large. 

Figure  1  shows  the  effect  of  a  single  point  Interference  with  a  small  circular  array.  It 
shows  that  if  the  interference  directions  differ  by  more  than  about  40  degrees  from  the  target 
direction  the  effect  on  F(n)  is  essentially  negligible.  It  must  be  borne  in  mind,  however,  that 
this  is  only  demonstrated  for  a  single  frequency  (5000  rad/ sec).  The  picture  looks  different  at 
other  frequencies,  and  the  integrated  effect  of  all  frequencies  therefore  has  the  effect  of  the 
loss  of  one  hydrophone  as  is  predicted  by  the  analysis  in  the  section  on  single  point  interference. 

Figures  2  and  3  are  similar  to  Fig.  1  except  that  the  interference  consists,  respectively,  of 
two  and  of  four  points,  separated  by  0.1  radian.  Since  the  interference  covers  a  large  azimuth 
segment,  the  effect  on  F(n)  covers  a  larger  angle;  however,  it  is  still  true  that  for  interference 
sources  at  angles  far  removed  from  zero  the  effect  on  F(n)  is  small.  A  similar  result  is  shown 
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Fig.  1.  Effect  of  single  point  interference 


Hn)  FOR  CIRCULAR  ARRAY  CONSISTING  OF  10  ELEMENTS  EQUALLY  SPACED  ON  A 


Fig.  2.  Effect  of  two  interference  points  spaced 
0.1  radians  apart  in  azimuth 
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Fig.  3.  Effect  of  four  interference  points  spaced 
0.1  radians  apart  in  azimuth 


Fig.  4.  Effect  of  two  interference  points  spaced 
90  degrees  apart  in  azimuth 
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F'(n)  Ft}}l  l  ist: Ail  Aim  AY  CONSISttNG  Of  10  ELEMENTS,  SPACED  2  FT  APART 


Fig.  5.  Effect  of  a  distributed  interference  source  approximated 
by  20  point  interferences  spaced  0.9  degree  apart.  Frequency  is 
4000  rad/sec. 


Fin)  FOR  LINEAR  ARRAY  CONSISTING  OF  10  ELEMENTS  SPACER  2  FT  APART 
FREQUENCY  IS  2000  HAD/SEC. 

INTERFERENCE  CONSISTS  OF  20  POINTS  LOCATED  AT  F  a  .  9°,  0,  .  1.8°, 
•■■‘,.17.1° 

K|M  ■  200 

ARRAT  IS  STEERED  ON  TARGET  AT  S  .  O0 


-80°  -W°  -<0°  4?  5  20°  80°  00°  80° 

INTERFERENCE  azimuth  ANCLE  F[ 


Fig.  6.  Effect  of  a  distributed  interference  source  approximated 
by  20  point  interferences  spaced  0.9  degree  apart.  Frequency  is 
2000  rad/sec. 
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In  Fig.  4  which  shows  the  effect  of  two  interferences  separated  by  a  large  angle  (30  degrees). 
The  last  two  figures  show  the  effect  of  a  strong  Interference  (K,  -  200)  distributed  over  a  rela¬ 
tively  large  angle  (17  degrees).  Again  the  effect  at  angles  far  removed  from  the  target  angle  is 
small,  but,  as  is  shown  in  Fig.  6,  the  relative  effect  is  quite  different  at  different  frequencies, 
as  has  already  been  pointed  out. 

The  computations  leading  to  the  results  shown  in  Figs.  1  through  6  were  quite  time  consum¬ 
ing,  with  computing  times  on  the  order  of  several  minutes  on  the  IBM  7094  for  the  cases  with 
large  numbers  of  interference  points.  For  this  reason  no  attempt  was  made  to  compute  the 
complete  detection  index,  since  this  would  have  required  summation  of  F(n)  over  a  large  num¬ 
ber  of  frequencies.  The  computer  results  therefore  still  do  not  conclusively  answer  the  ques¬ 
tion  of  how  serious  is  the  effect  of  large  distributed  Interferences.  The  indications  are,  how¬ 
ever,  that  the  results  of  the  previous  section  are  valid  under  considerably  wider  conditions 
than  those  assumed  there  to  produce  analytical  approximations.  In  fact,  it  appears  that  loss  of 
detectability  for  rather  widely  distributed  interference  Is  equivalent  to,  at  most,  a  few  hydro¬ 
phones. 


CONCLUSIONS 

The  major  difficulty  In  obtaining  general  estimates  of  the  effect  of  directional  noise  on  the 
detectability  in  an  array  processor  is  that  the  mathematical  manipulations  required  to  obtain 
answers  are  quite  complex.  Results  have  therefore  been  obtained  only  in  a  restricted  number 
of  simple  cases. 

The  general  tenor  of  these  results  is  that  if  the  ani&otropic-to-isotropic  noise  ratio  is 
small  the  effect  of  a  number  of  local  noise  sources  is  additive;  that  is,  the  loss  of  detectability 
resulting  from  two  noise  sources  of  equal  strength  is  twice  that  resulting  from  a  single  source. 
For  large  anisctropic-to- ambient  noise  ratio  the  effect  depends  on  whether  the  directional 
noise  sources  are  close  together  or  not.  For  a  single  point  source  it  has  been  shown  previously 
and  corroborated  here  that  the  loss  in  detectability  is  approximately  equivalent  to  the  loss  of 
one  hydrophone  from  the  array.  If  there  are  R  noise  sources,  widely  separated  from  each 
other  and  from  the  target  direction,  the  loss  is  approximately  equivalent  to  the  loss  of  R  hydro¬ 
phones,  provided  that  r  «  VST,  where  M  is  the  number  of  hydrophones. 

Point  noise  sources  that  are  close  together  affect  the  system  like  a  single  distributed  noise 
source,  and  the  indications  are  that  if  such  an  anisotropy  is  spread  over  a  relatively  small 
angle,  Us  effect  is  essentially  that  of  a  single  point  noise.  Unfortunately  this  has  not  been  con¬ 
clusively  demonstrated,  even  by  use  of  a  digital  computer,  and  only  a  rather  conservative  esti¬ 
mate  of  azimuth  angle  that  can  be  considered  to  be  "small”  has  been  obtained. 


.Appendix 

APPROXIMATE  INVERSION  OF  A  MATRIX  WHOSE  DIAGONAL  TERMS 
ARE  LARGE  RELATIVE  TO  THE  OFF-DIAGONAL  TERMS 


Let  the  n « n  nonsingular  matrix  A  be  given  by 

A  -  D+  B  ,  (A-l) 

where  D  is  diagonal  and  B  is  a  matrix  with  zero  diagonal  elements.  It  is  assumed  that  all  the 
non-diagonal  elements  of  B  are  of  about  the  same  order  of  magnitude,  and  that  the  elements  of 
D  are  of  about  K  times  that  magnitude,  with  K  »  i . 

The  inverse  of  A  is  given  by 

A' 1  =  (D+  B)'1  =  D-'(I  +  BD1)'1  =  D-‘  [l  -  BD' 1  +  (BD"  ')*  .  .  .  J  .  (A-2) 
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Since  the  elements  of  I)  are  of  order  K  relative  to  b,  the  elements  of  BD' 1  are  order  1  k  rela¬ 
tive  to  unity. 

It  can  be  shown  s  that  a  sufficient  condition  for  convergence  of  Eq.  (A-2)  is 


Ibjjl  <  1  .  i  . . . .  (A- 3) 


where  bSj  are  the  elements  of  bd-  1 .  Assuming  all  of  these  elements  to  be  of  about  the  same 
order  of  magnitude,  condition  (A- 3)  can  be  expressed  in  the  approximate  form 

nb0  l  ,  (A- 4) 

where  h0  is  a  representative  element  of  bd  This  element  is  of  order  1  K ;  therefore  con¬ 
vergence  requires 


n  K  <  1  .  (A- 5) 

The  convergence  will  clearly  be  more  rapid  if  this  inequality  is  sharper;  hence  one  can  ap¬ 
proximately  neglect  the  matrix  u  in  the  inversion  of  A  if  n/K  <<•  1. 
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